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Using a rotating flat layer heated from below as an example, we consider effects which lead to stabilizing an exponentially growing 
magnetic field in magnetostrophic convection in transition from the kinematic dynamo to the full non-linear dynamo. We present 
estimates of the energy redistribution over the spectrum and helicity quenching by the magnetic field. We also study the alignment of 
the velocity and magnetic fields. These regimes are similar to those in planetary dynamo simulations. 
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1 Introduction 

Many physical processes can be referred to threshold phenomena, when the increase of the governing pa- 
rameter leads to the appearance of an increasing solution. Such an example is thermal convection, when 
the growth of thermal instability starts at a certain critical value o f the Rayleigh number Ra'^'', which char- 
acterizes the amplitude of the heat sources ( Chandrasekhaiill96ll ). The same situation occurs if magnetic 
field B is generatied in a conductive medium: the increase of the convection inte nsity, paramet rized by the 
magnetic Reynolds number Rm, can lead to an exponentially growing solution ( Moffattlll978l ). After that 
magnetic field grows up to the moment when it starts to have an effect on the flow. As in many astrophys- 
ical objects Rm is very large, providing quite extended spectra of the fields, it is generally believed that 
this influence need not lead to a direct suppression of fluid motions. This statement is supported by the 
fact that in some cases transition from the non-magnetic to the magnetic state can be accompanied by 
the growth of Reynolds numbers. In other words, knowledge of Rm is not sufficient to answer the question 
whether the magnetic field will grow further or not. 

The most widespread point of view is that the magnetic field causes such a reconstruction of the flow 
that the generation of the mag netic field b ecomes less efficient. However, the visual control of the flow does 
not reveal an essential change ( Joneg|2000|), which can be due to the force- free nature of the magnetic field 

configurations: (V x B) x B 

One of the 



Lb/B^ <C 1, where Lb is the energy-carrying scale of the magnetic field, 
of 



1984 



explanations of the saturation mechanism is offered in (jBerger 

Brandenburg and SubramanianI |2005| ). where saturation is connected with the scale separation of 
the generated magnetic field due to conservation of magnetic helicity. In the present paper we study, 
on an example of the dynamo in a rapidly rotating flat layer heated from below, how such a magnetic 
energy redistribution over the spectra takes place in geostrophic systems. Our simulations demonstrate 
the occurrence of the magnetic a*^-effect which can suppress the total a-effect and lead to a saturated 
dynamo. 

The other point is the correlation of velocity and magnetic flelds (so-called cross helicity). It appears that 
a saturated velocity fleld can still lead to an exponentially growing magnetic fleld, provided that this new 
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artificial magnetic field does not contribute to the Lorentz force. It was first noticed by lCattaneo and TobiasI 
(J2009), see also thei r KITP's con f erence video presentation. This problem has now beco me the subject of 
various discussions ( Tilgneiil2008l . iTilgner and Brandenburg l2008l . ISchrinner et al.ll2009l ). We also present 
some results concerned with the magnetostrophic regimes close to those in geodynamo simulations that 
are also unstable for large Rm. 



2 Dynamo equations 

The geodynamo equations for an incompressible fluid (V • V = 0) in a layer of height L rotating with 
angular velocity Q in a Cartesian system of coordinates (x, y, z) in its traditional geodynamo dimensionless 
form can be expressed as follows: 

_ = V X (V X B) + q-^AB 



EPr" 



dV 
~dt 



+ (V-V)V 



-VP - I2 X V+ 



(1) 



dT 
'dt 



RaT zlz + (V X B) X B + EAV 



+ (V-V)(T + To) = AT. 



Velocity V, magnetic field B, pressure P and the typical diffusion time t are measured in units of k/L, 

K, is the thermal diffusivity, p is the density, /i the perme- 
is the Ekman number, v is the kinematic viscosity, r/ is the 



Y2fiK/ip;, pK /L and L /k, respectively, where k is the thermal diffusivity, p is the density, p the perme- 

K v 

ability, Pr = — is the Prandtl number, E = 

rrpT 

magnetic diffusivity, and q = kIt] is the Roberts number. Ra = is the modified Rayleigh number, 

a is the coefficient of volume expansion, ST is the unit of temperature, go is the gravitational acceleration, 
and To = 1 — z is the heating from below. The problem is close d with periodica.! boun dary conditions in the 
(x, y) plane. In the z-direction, we use simplified conditions ( Cattaneo et al.l 12003! ): for T and V: T = 

0, and the pseudo-vacuum boundary condition 



(heating from below), stress-free for V: V^ 



for B: B^ 



B„ 



M2 



dz dz 
at z = 0, 1. System ([1]) was solved using the pseudo-spectral Fortran-95 MPI 



code (JReshetnvak and Heidall2008l ) on cluster PC computers using grids N^, N = 64. 



3 Results of modelling 

3.1 General properties 

We consider two regimes that differ in Rayleigh and Roberts numbers: 

Rl: Ra = 4 • 10^, Pr = 1, E = 2 • 10~^ q = 10. 
R2: Ra = 1 • 10^, Pr = 1, E = 2 • 10-^ q = 3. 

To get the final saturated dynamo, we started from a pure convection state without the magnetic field. 

At time t = 0.35, we injected a magnetic field of small amplitude. After an intermediate kinematic regime 

^2 

with exponential growth of magnetic energy Em = , we arrive at a quasi-periodical state with kinetic. 



V 



2Ro' 



Ek = — , and magnetic, Em, energies of the same order of magnitude for Rl ( see Fig. [T]), and at the 
state with Ek > Em for R2 (see Fig. [2]), where a smaller value of q was used. 
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Figure 1. Regime Rl. Evolution of the mean over the volume fluctuations of the squared temperature T^, kinetic and magnetic energies 
Ek, Em, and kinetic helicity x^, Ix^l- 



Both regimes correspond to the geostrophic (magnetostrophic) state, see typical cyclonic structures of 
the temperature fluctuations T and kinetic energy Ek = — in Fig. [3]for R2. The diameter of the cyclones 

depends on the Ekman number as dc ~ LE^'^ ( Busselll97d ). In both cases the transition from the kinematic 
regime to the full dynamo was accompanied by an increase of the temperature fluctuations T^. 

On the other hand, the behaviour of the kinetic energy for the two regimes was different: for Rl the 
growth of the magnetic field leads to an increase in convection intensity {Ek becomes larger) and for R2 
Ek to a decrease. The corresponding Reynolds numbers are Re = 200 (kinematic dynamo regime) and 
Re = 420 (regime with saturation) for Rl and Re = 1500, 1100 for R2. For the observed geostrophic state 
the ratio of the nonlinear term in the Navier-Stokes equation to the Coriolis force, the Rossby number, 
Ro = EF is quite small: Ro = 4.3 10-^ 8.4 lO'^ for Rl and Ro = 3.0 lO'^, 2.210-2 for R2. One of the 
explanations of the growth of Re for regime Rl is that the magnetic field disturbs the regular geostrophic 
flow for Rl and helps the generation of the magnetic field. The total energy Ek + Em also increases. 
Regime R2 corresponds to the more disturbed state, and the kinetic energy is reduced to the amount of 
increase of the magnetic energy. This scenario usually takes place for large Reynolds numbers. Note also 
that, while for Rl the increase of T^ corresponds to the increase of the Archimedean work RaTT4, for R2 
the work decreases due to the chaotization of the fields. The behaviour of the kinetic helicity x = V • rot V 
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as a whole copies the evolution of the kinetic energy. The dispersion of x increases when the magnetic 
field reaches its saturated value. 



1.6 10 



f^ 8 10' 




Figure 2. The same as in Fig. [2] for regime R2. 



3.2 Spectra 

The maximum of the kinetic energy spectra for Rl corresponds to the horizontal scale of the cyclones 
kc = 1/dc, see Fig. Hl^l). The magnetic field slightly decreases this maximum and kinetic energy increases. 
The relative part of the kinetic energy on the large scales also increases. 

The growth of the magnetic field during transition from the kinematic dynamo state to saturation state 
is accompanied by an increase of the magnetic field on large scales Fig. [11^2). The first mode which reaches 
saturation is the mode with k ~ 10. The other modes still grow filling the spectra for Em- This behaviour 
is the same for both the regimes Rl, R2. The behaviour of the kinetic energy on the large scales is a little 
bit different for R2 ( Fig. [11^3)), where we observe a decrease of Ek, which corresponds to the breakup of 
the horizontal rolls by the magnetic field. The maximum for Ek at k = 1 then disappears. 
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The inhomogeneous growth of the magnetic field for different k is quite important for understanding 
the saturation mechanism of the magnetic field. The growth of the magnetic field in t he kinemat i c state 
takes place at convective times r^ ~ {kVk)~^, which decreases with k. According to iKazantsevI ( 19681 ) 
spectrum of the magnetic field is Em ~ k^''^ for the non-rotating turbulence, and the maximum of the 
magnetic field is then close to the dissipative scales. In our regime, the maximum of Ex at k = kc is more 
important for the magnetic field distribution, and the modes with k ^ kc reach saturation level at first. 
We argue that this regime i s closely connected with the occurrence of the coherent structures discussed in 
(iTobias and Cattaneolblx)^ ). 




-0^07 



Figure 3. Distribution of temperature fluctuations T and kinetic energy Ek for regime R2. 

The observed redistribution of the magnetic field o ver the scale s is closely related to the mean over the 



volume magnetic helicityy-'^ = A • B for Rm — >■ oo (;Bergeii Il984l ) , where A is the vector potential of the 
magnetic field B = rotAJ^. T he alternating -sign quantity x"^ is a measure of the linkage of the magnetic 



field lines with one another (|Moffattlll978l ). After so me algebra with the induction equation, one arriv es 
at the equation for the mean over the volume fields ( Bergeiill984l . iBrandenburg and Subramanianll2005l ) 



^ A • B = -R^i J • B + n. 



(2) 



where J = rotB is the current, and 11 is the flux of x through the boundary. 11 = for the fully periodical 
boundary conditions as well as for the super conductive boundaries . The other scalar product x = J • B 
is the so-called current helicity. 

Then, after time t ~ Km ^ 1, one has saturation regime -— - = 0, with 



J B 



-jb. 



A B 



-a • b, 



(3) 



where the decomposition of field F into the mean and fiuctuating parts was used: F = F + f 
( Brandenburg and Subramanianll2005l ) . In other words, after the kinematic regime no changes of the mag- 
netic field would change the total magnetic and current helicities. Any local change is possible only due 



^ Large Rm is typical for many astrophysical bodies which possess their own magnetic field. 
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Figure 4. Spectra of the kinetic and magnetic energies for regime Rl (1-2) and R2 (3-4). The solid line correspond to the saturated 
dynamo, squares to the convection without magnetic field, and crosses to the kinematic dynamo regime. 

to the redistribution of x^i X^ over the scales. Of course, this approach becomes more complicated, if 
the mean helicities change sign in space, and the ideas on scale decomposition are applicable to the space 
domain with the same sign of helicity. 

The application of the pseudo-vacuum boundary conditions, which provides 11 ^ at the bound- 
aries, is more t ricky. These conditions l eads to the increase of the energy of the magnetic field 
at large scale s ("Brandenbure: a nd SandinI 12004 ) and break the catastrophic quenching predicted by 
Vainshtein and Cattaneo ( 1992) a nd observ ed for the fully periodic boundary conditions, see (JHughes et al. 



199d. l Brandenb ur. and Subramanll^ Ei^^. There is also an indica tion of catastrophic gu enchin. for the 



mean- field dynamo models with periodic boundary conditions, see ( Cattaneo and Hugheslll996l ). Accord- 



would be saturated at Bq 



R^^5^, where 6^ 



~ V 



ingly to this scenario the large-scale magnetic field Bq 
are the small-scale field energies. 

In our simulations with pseudo-vacuum boundary conditions, we still observe some decrease of the 
magnetic field intensity, see Fig. [THSJ the ratio of the magnetic to kinetic energies decreases from Rl to R2, 
which should be explained by a decrease of the Roberts number q rather than by catastrophic quenching. 
In Fig. Owe also observe that the magnetic helicity on the large scales increases sufficiently after transition 
to the saturated regime . 

Let us recall that, according to the mean field dynamo theory, there is connection between the ki- 
netic and current he l icities and hydrodyna mic a - and magnetic a -effects for short correlation times 
(IPouquet et al.lll97()l . IZeldovich etfaPligS^I ) 



a 



U 



-TVT7J/3, a-^ = Tj • b/3, (4) 

where u = rotv is the vorticity and r is a correlation time. In practice ( Zeldovich et al.lll983l ) these formulas 
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Figure 5. Spectra of the current helicity for regime Rl. The sohd line corresponds to the saturated dynamo, circles to the kinematic 
dynamo regime. 
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Figure 6. Z-profiles of the kinetic x^ (solid line) and current x^ (squares) helicities for regimes Rl (1) and R2 (2). 

are applicable when the typical time of the large-scale magnetic field growth tl is larger than the turnover 
kinetic time ti = l/v{l), where v{l) is a velocity on scale I. A rough estimate of tl, using jump of the 
magnetic energy at the kinematic regime Rl, (see Fig. [1]), yields tl ~ 10^^, which is already smaller than 
Tt ~ 1/ (Refcc) ~ 10~^ (here we supposed that the kinetic energy is concentrated in the vicinity of kc)- 
Taking into account that the mean magnetic field grows slower than the small-scale field, we find that tl 
is even much larger than the above estimate (a similar situation for regime R2 takes place). 
The total a-effect is then 



a 



a-^ + a^. 



(5) 



If the signs of the helicities are the same, the total a-effect is reduced (a — )■ —a") and the magnetic field 



stops growing. The latter is well observed in our simulations. Fig. [6HZ1 where x > ^^r z < 0.5 and 

X < for z > 0.5. All the three helicities x have the same signs, however, x is more rugged 

due to the contribution of the small-scale fields. 



3.3 Energy fluxes in the wave space 

Energy redistribution over the spectra is closely related to the fluxes in wave space. In spite of the fact 
that the final saturated state is quasi-stationary, there are still flu xes in wave space. Thi s happens because 
the scales where energy is generated and dissipated are different ( Rose and Sulemlll978l ). For the simplest 
cases, such as 3D Kolmogorov's turbulence, k inetic energy goes from the scale of the force to the the 
dissipative scale. Here we repeat the results of (JReshetnyak and Hejdall2008l . lHejda and Reshetnyakll2009l ) 
and discuss the difference between energy fluxes of kinematic dynamo regime and the saturated regime. 
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Figure 7. Z-profilos of the magnetic x'^ ^nd cross- x helicities for regime Rl (1-2) and R2 (3-4). 



d 
Consider the kinetic energy flux in wave space through wave number k: Txik) = —-^ < (V x rotV) • 

^ CO 

V^ >, where V<(r) = ^ Vp e*^"" is a low frequency counterpart, and f Txik) dk = 0. 

|P|<fc fc=0 

Using relation rot (V x B) = — (V • V) B + (B • V) V, the magnetic energy flux TM{k) can be decom- 

oo 

posed into advective and generating parts Tm = Tn — Tl, J TN{k) dk = 0. The latter is equal to minus 

fc=o 
the work of the Lorentz force. 

The fluxes of kinetic energy Tk are presented in Fig. EI^l), [9{1). For k ^^ kc the inverse cascade of the 
kinetic energy is observed: the cyclones are the sources of the energy. Here energy is distributed to the 
large scales {k < kc, Tk > 0, inverse cascade), as well as to the small scales {k > kc, Tk < 0, direct 
cascade) where it dissipates. The magnetic field causes some blurring of the maxima and shift of Tk to 
the large-scale region. 

In contrast to Tk, Tm includes not only an advective term, but also a generating term. It appears that 
the integral of 7m over all k is positive. Moreover, Tm is positive for any k, Fig. [5[2),[U|^2), i.e. the magnetic 
field is generated on small scales. Position of the maximum of Tm is close to the maximum in the spectrum 
of Em- The form of Tm is the same for the small magnetic field during kinematic regime and for the 
saturated mode. Let us consider where the magnetic energy comes from: Is the energy transfered from 
other scales to some certain k, or is it generated localy? Generating flux —71 is shown in Figs [11^3), [91^3). 
Its maximum is close to the minimum of Tk (Rl); i-e the kinetic energy is transformed to the energy of 
the magnetic field. Except for a small region at high k, there is always an inverse cascade of magnetic 
energy —71 > 0. Note that the amplitudes of the maximum of —71 and Tm are sufficiently different. That 
is because —71 and Tn are anti-correlated, see Fig. [8]J4), [9^4). This means that advective term Tn < 
transports the major part of the energy to the dissipative scale. 

Transition to the saturated state for Rl does not change form of the flux's curve substantially. It only 
increases the amplitudes of —71 and Tn, thus increasing the synchronization of the fluxes —71 and Tn- The 
maxima of —71 and Tn shift to the large scales for the full dynamo regime. For the R2-regime saturation 
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Figure 8. Fluxes in the wave space of the kinetic Tk and magnetic Tm energies, Lorentz work Tl and the convective term in the 
induction equation Tpf, regime Rl. The solid line corresponds to the saturated dynamo and crosses to the kinematic dynamo regime. 



changes the flux of the magnetic energy more efficiently. For the kinematic regime, the flux was more or 
less homogeneous from the small to large scales. For the saturated state, cyclones at kc conected with 
decrease of Em/Ek provide the main contribution to — 7l- The relative strength of the total flux of the 
magnetic energy due to the —71 term to the main scale remains at the same level. This flux can be related 
to the a-effect. 



4 Alignment of the fields 

In discussing z-profiles of the fields, we have not yet taken into account the importance of the spatial- 
temporal correlation of fields V and B. In this connection an interesting question arises: Can an already 
quenched velocity field in the dynamo model generate an exponentially growing magnetic field or not 
( Cattaneo and Tobias! l2009l . iTilgner and Brandenburg) l2008l )? The only one difference between this new 
passive magnetic field and the original one is that the new field does not contribute to the Lorentz force. 
It appears that the answer depends on the spectrum of the magnetic field and its time behaviour. Usually , 
if one has only one excited magnetic field mode for the saturated regime (ITilgner and Brandenburdl2008l) 
the new, passive magnetic fi eld is stable. This statement is supported by simulations in the sphere (ITilgner 
20081 . ISchrinner et al.ll2009l ). where the dipole mode dominates for low Rossby number Rq . For the multi- 
mode regime (which is the case for the large Rayleigh number) the situation is different: the new magnetic 
field B grows exponentially ( Cattaneo and Tobiasll2009l ) . There is an indication that the threshold for such 
a different behaviour of B in time in the sphere takes place at Rq = 0.12 ( Schrinner et al.ll2009l ). Here we 
consider what happens in the fiat layer dynamo for regimes Rl, R2 with quite extended spectra, adding a 
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Figure 9. The same as in Fig. [8] for regime R2. 

new induction equation for B to the system ([1]): 

-— = rot (V X B) + q^^AB 



EPr~ 



dT 



dV 
'dt 



+ (V-V)V 



-VP - Iz X V+ 



RaT zlz + (V X B) X B + EAV 



(6) 



— + (V-V)(T + To) = AT 



rot ( V X B ) + q^^AB 



dB 

'dt 

For both the regimes we observe distinct behaviours of B^ and B^, see Fig. [lOl Although the same 
velocity field V was used in both the equations for B and B , the new artificial magnetic field starts to 
grow exponentially. We are ready to conclude that temporal synchronization of the fields in space and time 
is crucial for stabilization. Field B is certainly less synchronized with velocity field V because the Lorentz 
force based on B has been omitted. This idea is supported by our simulations of the autocorrelation 
functions of the magnetic and velocity fields Corr{Vi, Bi), Corr{Vi, Bi) calculated over half of the volume 
z < 0.5, see Fig. [TT| [T2l Evidently for both the regimes Rl, R2 the correlation for B is less than for 
field B. However, the correlation is quite small for both the regimes due to the stochastic nature of 
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Figure 10. Evolution of the magnetic energies _B^/2 (solid line) and B^ /2 (circles) for regimes Rl (1) and R2(2). 
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Figure 11. Evolution of the autocorrelation functions for regime Rl. The thin line corresponds to the original magnetic field B and the 
line with circles to B. 



the small scale fields and does not exceed a few percent, reducing with increasing Rm. The reduction of 
correlation for all the components corresponds to the reduction of alignment of fields V and B, when the 
Lorentz force is omitted. The real magnetic field B has strong er alignment. We adopt the explanation 
given bv lCattaneo and Tobiad ( 2003 ) . iTilgner and Brandenburg] (|2008: ) that the instability of B is due to 
the difference in the stability criteria for the full dynamo equation with quenched Vq, and single induction 
equation with given Vq. 
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Figure 12. The same as in Fig. llll for regime R2. 



5 Discussion 



The transition from kinematic to saturated dynamo regime in cyclonic convection is accompanied by the 
reconstruction of the flow as well as of the magnetic field. In general, the change of the kinetic energy is 
not crucial: it can even increase during the transition. More important is the reconstruction of the flow 
patterns. The growing of the magnetic field B from the quasi-stationary convective state with non-zero 
kinetic helicity x {^) is defined by the leading eigen solution. This first mode grows up to the level when 
the maximum of the spectra at k r^ kc reaches the saturated level close to the equipartition value on 
this scale. The growth of field B then stops at kc- The small-scale field produces the magnetic a -effect 
which suppresses the total a-effect ([5]). According to ([3]), the transition to the saturated regime requires 
an increase of the large-scale magnetic field, which takes place at diffusion time t ~ Rm- Some change 
of the magnetic energy flux in the wave-space is also observed. The long-term fitting to the saturated 
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regime is also predicted by the dynamical models of a-quenching ( Kleeorin et alJll995l ). During this time 



the large-scale magnetic field still grows to the limit defined by the restriction on the magnetic helicity 
conservation. The Lorentz force, which provides the correlation of the velocity and magnetic fields V and 
B and their alignment, plays a crucial role. 
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